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Abstract 

In order to prove irrationality of v2 by using only decimal expansions 
(and not fractions), we develop in detail a model of real numbers based 
on infinite decimals and arithmetic operations with them. 



1 Introduction 

Irrationality of V2 via decimals? Many proofs of irrationality of the number 
y/2 were published and collected, e.g., Beigel [4], Bogomolny [5] (at least 19 
on-line proofs), Flannery [IB], Gardner [TH] (lists 18 references), Harris [23] (13 
proofs) , Miller and Montague [32] , Myerson [33J , Subbarao [32] , and Waterhouse 
[44] , and many textbooks on algebra, mathematical analysis and number theory 
present such proof: Allouche and Shallit [U Theorem 2.2.1], Apostol [21 Theorem 
1.10], Hardy and Wright [23j Theorem 43], Jarm'k [26j str. 16], Pugh [34] p. 
11], Tao [33J Proposition 4.4.4], Zorich [33 p. 51], to name a few. These proofs 
start invariably from the hypothetical fractional representation \f2 = ^ with 
a, b € Z and derive a contradiction. But could one prove that 

V2 = 1.414213562373095048801688724209698 . . . 

is irrational in a purely numerical way by showing that this decimal expansion 
is not ultimately periodic, without invoking the fractional model of rational 
numbers? Such proofs are hard to come by in the literature. We know of only 
two imperfect examples and reproduce one of them below. (The other is proof 
no. 12 in [5], due to A. Cooper, which is less convincing as it fails to deal, 
without help of fractions, with the case 2 below.) 

A proof of irrationality of \/2 by decimals. Lindstrom [29] found the 
following proof of irrationality of y/2 using decimal expansions. It brings to 
contradiction the assumption that \[2 is a terminating or repeating decimal. 
We abridge it somewhat but preserve its style and notation. 
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Case 1. If V2 is a terminating decimal, say \/2 — l.abc with decimal 
digits a, b, c and c is nonzero, multiplying by 1000 and squaring, we 
get {labc is in the decimal notation) 

2000000 = (1000V2) 2 = {labc) 2 . 

This cannot hold as the units digit of 2000000 is zero but the units 
digit of the right-hand side is nonzero, the units digit of c 2 . Other 
terminating decimals are treated in a similar way. 

Case 2. If V2 is a repeating decimal, say \[2 = l.abc = l.abcabc . . . 
with some decimal digits a, b, c, then {abc is in the decimal notation) 
\/2 = 1 + abc/ 10 3 + a6c/10 6 + . . . . Summing the geometric series, 
rearranging and squaring, we get 

2 • 999 2 = (V2(10 3 - l)) 2 = ((10 3 - 1) + abc) 2 . 

This cannot hold either as the units digit on the left-hand side af- 
ter multiplying is 2 but the units digit on the right-hand side after 
squaring is that of a square, one of 0, 1,4,5,6,9. Other repeating 
decimals are treated in a similar way. 

Hence V2 is irrational. 

Is anything wrong with the proof? It is not and it is. In the liberal 
view that allows mixing infinite decimals with fractions, the above proof is 
fine. Indeed, to be fair, it appears that its author did not set as his goal to 
use exclusively decimal expansions. But in the more restricted view when one 
accepts only arguments using decimal expansions, the above proof suffers from 
two problems. First, in summing 1 + abc/10 3 + abc/10 e + ... it uses fractions. 
Second, and this is a more serious problem, it uses the familiar but in our context 
problematic simplification (V2) 2 = 2. When real numbers are infinite decimals, 
the correct simplification is (V2) 2 = {1.999 2}, i.e., one has to consider 
two possibilities, (\/2) 2 = 1.999... and (V2) 2 = 2. Indeed, (1.4142135 .. .) 2 , 
calculated via truncations as 1.4 2 , 1.41 2 , 1.414 2 and so on, approaches 1.999 . . . 
and not 2. The above proof does not consider the possibility (\/2) 2 = 1.999 . . . 
at all and therefore is incomplete. 

But wait, may object the attentive reader, isn't it so that we identify 1.999 . . . 
with 2 and thus it suffices to work just with the representative 2? Therefore, 
isn't the argument in case 2 looking only at (V%) 2 = 2 sufficient after all? It 
is and it is not. It is sufficient because when the model of K based on decimal 
expansions is worked out in detail (in Section 2), it turns out that it really 
suffices to look, in a sense, only at the possibility (V2) 2 = 2 (see the proofs 
of Propositions 13.11 and I3.3[) . At the same time it is not sufficient because if 
such model is not specified or at least referred to, which, alas, is the rule in 
discussions on this topic, all arguments about what can or cannot be justified 
by identifications like 1.999 . . . = 2 are necessarily scholastic and lack substance 
(cf. various "proofs" that 0.999 . . . = 1 at [35]). 
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What is in this text. Our original motivation was to understand if the proof 
of Lindstrom qualifies as a "numerical" proof of irrationality of \f2 using only 
decimal expansions and if it does not, how would such a proof look like. As to 
Lindstrom's proof, we summarized our opinion in the two preceding paragraphs. 
In Section 3 we present our proofs of irrationality of y2 using only decimal 
expansions. They require a model of R based on decimal expansions, called a 
decimal model o/R for short, which we develop in detail in Section 2. 

Such model is of course not a complete novum because besides the two 
well known 1872 models of R, the model of Meray, Cantor and Heine using 
fundamental sequences of rationals and Dcdckind's model based on cuts on 
rational numbers, there is, allegedly, a third model based on decimal expansions, 
developed by Weierstrass and Stolz. Unfortunately, unlike for the first two 
classical models, we could not find in the literature any really detailed and 
satisfactory presentation of the Weierstrass-Stolz model or its modern version, 
only outlines and sketches; we give references and more comments at the end of 
Section 2 where we present our stab at this model. 

Troubles with formal limits. So how to add and multiply infinite decimals? 
Fowler [17] challenges the reader to calculate the product 

1.222222... x 0.818181... = ? 

A natural idea is to take finite truncations of the factors, form the sequence of 
their products, 

1.2 x 0.8 = 0.96, 1.22 x 0.81 = 0.9882, 1.222 x 0.818 = 0.999596, 

and take its formal limit. Each digit after the decimal point stabilizes eventually 
at 9 and each digit before it is always 0. Hence, indisputably, 

1.222222 ... x 0.818181 . . . = 0.999999 .... 

In this way, using formal convergence of digits, one can define multiplication 
and addition of infinite decimals. 

This often suggested approach brings certain troubles. It takes an effort to 
show that subtraction is well defined because the sequence of differences of the 
truncations is in general not eventually monotone. But the main trouble is that 
formal limits do not commute with arithmetic operations and consequently these 
operations on infinite decimals lose their convenient properties. The following 
examples show that addition ceases to be associative and that the distributive 
law fails: 

(-0.999... + 1)+ 0.999... = 0.999... but - 0.999 ... + (1 + 0.999 ...) = 1 
and 

(10- 1) • 0.111 ... = 0.999. . . but 10 • 0.111 1 • 0.111 ... = 1. 
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Certainly, when identifications like 0.999 . . . = 1 are applied these irregularities 
go away (as one expects they must) but again it takes certain effort to prove it 
rigorously (without handwaving). 

We compare commutativity and associativity of arithmetic operations on 
infinite decimals. Since these operations on the finite truncations are clearly 
commutative, applying formal limits we get that addition and multiplication of 
infinite decimals is commutative. At first one might think that, similarly, so 
it works for associativity (cf. Gowers |22| ) but, with the previous example, we 
know very well that this argument is erroneous. Addition of finite decimal trun- 
cations is of course associative but formal limits do not transfer this to infinite 
decimals. The point is that (unlike commutativity) associativity of an operation 
with infinite decimals amounts to the exchangeability of two formal limits, which 
is a nontrivial result to be proven, and consequently is not at all automatic from 
the associativity for finite truncations. For multiplication this exchangeability 
occurs and multiplication of infinite decimals is associative. For addition it fails 
and addition is not associative; the same happens for the distributive law. (The 
reason is that for multiplication one has always monotonicity of partial products 
of truncations, which is in general not the case for subtraction.) 

Let us also illustrate the difference between usual metric limit and formal 
limit of decimal expansions by a curious paragraph from Courant and Robbins 
[TUl p. 293]. After explaining the standard notion of the limit a of a real 
sequence ai, a^, ■ . ■ , they write 

If the members of the sequence oi, 02, a 3 , • • • are expressed as infinite 
decimals, then the statement lima„ = a simply means that for any 
positive integer m the first m digits of a n coincide with the first 
m digits of the infinite decimal expansion of the fixed number a, 
provided that n is chosen sufficiently large, say greater than or equal 
to some value N (depending on m). This merely corresponds to 
choices of e in the form 10 _m . 

So, since the first m digits of a n are for large n equal to the first m digits of 
the decimal expansion of a, should one conclude that each digit of a n eventually 
stabilizes and the decimal expansions of a n formally converge? It is not exactly 
what the paragraph says (though it may seem to follow immediately from it) and 
it is not actually true, as shows the sequence a\ = 1.1, a>2 = 0.9, 03 = 1.01, 0,4 = 
0.99,a5 = 1.001,06 = 0.999, . . . that metrically converges to a — 1 but has no 
formal limit (replacements of 0.999 by 0.998999 ... or 1.001 by 1.000999 . . . etc. 
change nothing on this). The correct reading is that the first m digits of a n 
are for large n equal to the first m digits of a decimal expansion of a. In our 
example, one needs to switch constantly between the two expansions 0.999 . . . 
and 1.000 ... of a = 1, the former being used for ai n and the latter for am-\- Is 
this what the authors meant? (No further explanation or example are given.) 

Our approach. We had tried formal limits as well but then abandoned them 
in favor of an alternative approach. Instead of developing the cumbersome 
arithmetics of infinite decimals and fixing it at the end by identifications like 
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0.999 . . . = 1, we work with them from the very beginning and define arithmetic 
operations with infinite decimals by means of a limit weaker than the formal 
limit as inherently multivalued (i.e., sometimes bivalued) operations. Our an- 
swer to Fowler's question is that 

1.222222 ... x 0.818181 . . . = {0.999999 . . . , 1} 

and our definition of multiplication (and addition) is such that both possibilities 
for the result are on completely equal footing. Similarly we define, for example, 
that 

{0.1999 . . . , 0.2} + {-0.5, -0.4999 . . .} = {-0.3, -0.2999 . . .}. 

This is a less ambitious approach than formal limits as we do not care which 
of the two possibilities of the result is the "correct" one. Its advantage is that 
all required properties of arithmetic operations with real numbers can be estab- 
lished in a straightforward and natural way and thus, finally, a sound numerical 
proof of the irrationality of y/2 can be formulated. However, at the end of Sec- 
tion 2 we return to formal limits and prove that they do provide well defined 
arithmetic operations with infinite decimals. 

2 A decimal model of R 

We start from the ordered ring of integers Z = (Z, +,-,<), which we assume 
to be given, and build from it the complete ordered field of real numbers M = 
(R, +, •, <) and, on the way, the field of rational numbers Q = (Q, +, •). By a 
ring we always mean a commutative ring with 1. We use notation N = {1, 2, . . .} 
(natural numbers) and Z = {. . . , —2,-1, 0, 1,2,.. .} (the integers). 

Decimals and real numbers. A signed decimal d is a pair of the sign + 
or — and an infinite string of decadic digits which are indexed by the integers 
k, k — 1, k — 2, . . . and the first of which is nonzero: 

d = ±cikCik-iak-2 ■ ■ ■ 

where k £ Z, a,- L e {0, 1, 2, 3, 4, 5, 6, 7, 8, 9} C Z and a k ^ 0. The zero decimal 
has no sign, k = 0, and has only zero digits: 

= aoa_ia_2 . . . = 000 .... 

The set of all decimals, signed and zero, is denoted B. Decimals with the + 
sign are positive and those with the — sign negative. For d € B and i e Z we 
write di for the zth digit at of d. If at does not exist (i > k for signed decimals 
or i > for the zero decimal) we set di = 0. The equality of two decimals d 
and e means that they have equal signs and di — ej for every i e Z, or both are 
zero. Note that we regard digits in decimals as integers (not just labels) and 
add, multiply and compare them as elements of Z. The decimal differing from 
a decimal d only in sign is denoted — d (for the zero decimal we set —0 = 0). 
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When writing down decimals we use usual conventions concerning the decimal 
point and omitting the + sign and the trailing zeros. 

We define a linear order (D, <) on decimals. It is the lexicographic ordering, 
adjusted for signs. More precisely, for d, e e D, d ^ e, we set d < e if and only 
if d and e have the same + (respectively — ) sign and there is an m £ Z such 
that di = e, for i > m but d m < e m (respectively d m > e m ), or d has — sign 
and e = or e has + sign, or d = and e has + sign. 

This linear ordering is not dense, for some pairs of decimals d < e there is 
no decimal / with d < f < e. It happens exactly when d and e have the same 
+ (respectively — ) sign and there is an m £ Z such that di — e, for i > m, 
e m — d m = 1 (respectively e m — d m = —1), and for i < m one has di = 9 and 
e j = (respectively d t — and a — 9) . We call such pairs d < e jumps. Jumps 
are clearly disjoint. Examples of jumps are 

0.999... <1 or - 17.341 < -17.340999... . 

We write T9 for the decimals ending with infinitely many 9's and To for the dec- 
imals ending with infinitely many 0's. We call the latter terminating decimals. 
Note that 6 To. Every jump consists of one element from To and one element 
from Tg and every element of To U Tg\{0} appears in exactly one jump. 

Let ~ be the equivalence relation on D identifying elements in jumps, i.e., 
c?~eiffd = eorifc?<eore<disa jump. We define the set of real numbers 
as the set of equivalence classes. 

Definition 2.1 The set of real numbers R is 

E = E>/~ = {{0}}U {{d} I de D\(T UT 9 )}U {{c,d} \c<d is a jump}. 

For IcDwe let [X] denote the set of elements of R intersecting X; for d E D 
we write [d] (and not [{d}}) for the equivalence class of d. For a, (3 € R we set 
a < (3 iff d < e for some decimals d £ a and e £ (3. It is clear that (R, <) is a 
dense linear ordering. 

We shall prove the following two theorems. 

Theorem 2.2 The structure R = (R, +,-,<), with the operations + and ■ on 
the real numbers still to be defined, is a complete ordered field. 

We let P denote the set of ultimately periodic decimals: d £ P iff there are 
m £ Z and p £ N such that di — di- p for every i < m. We have To, Tg C P. 

Theorem 2.3 The prime field of the field R = (R, +, •), isomorphic to Q, is 
formed exactly by the ultimately periodic real numbers [P] . 

Before we introduce arithmetics of real numbers, we dispose with the com- 
pleteness in Theorem [221 and show that in (R, <) every nonempty set bounded 
from above has supremum, i.e., the least upper bound. It is clear from how 
(R, <) arises from (D, <) that it suffices to prove this for decimals. 
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Proposition 2.4 In (D, <) every nonempty set bounded from above has the 
least upper bound. 

Proof. Let D C D, D ^ 0, be bounded from above. We may assume that all 
elements in D are positive; other cases easily reduce to this. Since D has an 
upper bound, there is a k G Z such that d, = for every d € D and i > k 
but dfe > for some d E D. We set = max (4, taken over all d G -D. Hence 
Cfe > 0. We set c^-i = max4-i, taken over all d € D with dfc = c k - We set 
Cfc_2 = max<ifc_2, taken over all d G Z? with d^ = Ch,dk-i = Ck-i- Continuing 
this way, we define a decimal c e D, with + sign. It is immediate that c is the 
least upper bound of D. □ 

Similarly for infima, i.e., largest lower bounds. 

Arithmetic operations with real numbers. We get them as extensions of 
the unproblematic arithmetic operations with terminating decimals, which we 
begin with. Let Z be the set of all decimals d satisfying di = for every i < 0. 
The mapping (decadic notation) 0^0 and 

±a fe a fe _i . . . a 000 . . . i-» ±(a fe 10 fc + a fe _ilO fe_1 H h a 10°) 

is a 1-1 correspondence between Z and Z. We get an ordered ring (Z , +, ■, <), 
isomorphic to the ring of integers. For k e Z we let 10 fc denote the positive 
terminating decimal with 

(10 fe ) fc = 1 and (10 fc ), = for i ^ k. 

For d G D and k e Z we denote by 10 fe d the decimal obtained from d by keeping 
its sign and shifting its digits by k places: 

(10 fc d) 4 = d^ k . 

Note that 10 fc (lCT fe d) = d for every k € Z and deB and that the elements d 
and 10 fe d of Zo are mapped by decadic notation to the elements a and I0 k a of 
Z. For every d e T there is an N e N such that 10™d e Z for every n > N. 
To add and multiply two terminating decimals c and d, we take some fc e Z 
such that 10 fc c and 10 fe d lie in Z and define 

c + d= l(T fe (10 fe c+ 10 fe d) and c • d = l(T 2fe (10 fe c • 10 fc d) 

where the operations on the right sides are in Z . It follows that the results do 
not depend on k. Also, for c, d e T we have c < d iff the terminating decimal 
d — c has + sign (i.e., the order on Z transferred from Z coincides with the 
already defined lexicographic ordering). 

Proposition 2.5 (To, +,-,<) is an ordered ring. 
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Proof. The elements and 1 (= 1.000 . . .) are neutral to addition and multi- 
plication and any d £ To has additive inverse —d. The required properties of 
+, • and < (commutativity, associativity, distributivity and monotonicity) follow 
from the fact that they hold in Zo. For example, if c(d + e) ^ cd + ce for some 
c, d, e £ T , we take a k £ N so that all 10 fc c, 10 k d and 10 fe e lie in Z and obtain a 
refutation of the distributive law in Z : 10 2fe (c(d + e)) = (10 fe c)((10 fe d) + (10 fc e)) 
differs from 10 2fe (cd + ce) = (10 fe c)(10 fc d) + (10 fe c)(10 fc e). Similarly for other 
properties. □ 

Two defects of this ring, non-completeness and the lack of division, will be fixed 
by extension of the operations to K. 

For n £ N and d £ D, the n-truncation d\n £ T has the same sign as d and 
digits 

(d\n)i — di for i > —n and (d\n)i = for i < —n. 

If (d\n)i = for every i £ Z, we omit the sign and set d\n to be the zero 
decimal. Any truncation of the zero decimal is the zero decimal. For a sequence 
of terminating decimals c(n), n = 1, 2, . . . , we write 

c(n) -» 

if for every k £ N there is an N £ N such that c(n)\k = whenever n > N. 
It is clear that c(n) — > and d(n) — » implies c(n) + d(n) — » 0. Similarly, 
c(n) — c'(n) — > and d(n) — d'(n) — + for four bounded sequences of terminating 
decimals implies c(n)d(n) — c'(n)d'(n) — ► 0; it follows from the rearrangement 

c(n)d(n) - c'(n)d'(n) = (c(n) - c'(n))d(n) + c'(n)(d(n) - d'{n)). 

The next key result gives an alternative arithmetic characterization of jumps. 

Proposition 2.6 Two decimals c and d are equivalent, c <~ d, if and only if 
c\n — d\n — > 0. 

Proof. Suppose that c ^ d. If c = d then for every n even c|n — d|n = 0. If c < d 
or d < c is a jump then c|n — d|n = ±10 _n for large n and so c|n — d\n — > 0. 
Suppose that c ^ d. We assume that d < c is not a jump (the other case c < d 
is similar). So d < e < c for a decimal e. We have that e\n — d\n > 10~" for 
every large n and the same holds for c\n — e\n. Thus there is an m £ N such 
that c|m — d|m > 2 • 10~ m . But then for every j > m we have that 

c|j - d|j > 2 • 10~ m - 9 • lO-" 1 " 1 - 9 • 10-" 1 ^ 2 9 • 10" j > 10~ m , 

which means that (c\j — d\j)\m ^ 0. Thus c|n — d\n ■/* 0. □ 

Decimals in a jump are distinct but infinitesimally close; this phenomenon occurs 
here without any nonstandard analysis. 

We say that c <G D is a hybrid limit of a sequence c(n) of terminating decimals, 
written hlim„ c(n) — c, if 

c(n) — c\n — > 0. 
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In other words, hlim n c(n) = c means that for every k G N there is an N G N 
such that \c(n) — c\n\ < 10~ fc whenever n > N. Hybrid limits in general arc 
not unique: each of the three sequences c(n) — 1, c(n) = 0.99 . . . 9 = 1 — 10~ n , 
and c(l) = 1, c(2) = 0.9, c(3) = 1, c(4) = 0.99, c(5) = 1, ... has hlim„ c(n) = 1 
and hlim n c(n) = 0.999 . . .; in the formal sense the first sequence converges to 
1, the second to 0.999 . . . and the third is not formally convergent. To make 
it unique, we will regard hlim„ c(n) as the set of all its (at most two) values. 
By Proposition 12.61 hlim„ c(n) = c and hlim„ c(n) = d implies c ~ d, and 
hlim n c(n) = c and c ~ d implies hlim„ c{n) = d. So hlim„ c(n), if it exists, is 
a unique element of R: 

hlim„ c(n) = or hlim„ c(n) = {c. d} = [c] = [d] G R 

(here c = d or c ~ d). Note that if c(n), d(n) € To with c(n) — d(n) — > then 
hlim„ c(n) = hlim ra d(n). We are ready to define addition and multiplication of 
real numbers. 

Definition 2.7 For a, (3 G R we set a + (3 = hlim„ (c|n + d\n) and a/3 = 
hlim n (c|n • d\n), where c G a and d G are arbitrary decimals. 

By the previous remarks, the results are independent of the selection of c and 
d, and the resulting limits, if they exist, are real numbers. This completes the 
statement of Theorem 12.21 and it remains to prove that the defined arithmetic 
operations are always defined and have required properties. 

Properties of formal and hybrid limits. To accomplish it we first derive a 
few properties of hlim n . It helps to use also formal limits. We say that c G D is a 
formal limit of a sequence of terminating decimals c(n) , written flim„ c(n) = c, 
if for every k G N there is an N G N such that for every n > N the decimals c 
and c(n) have same sign or c is zero, and 

c(n)\k = c\k. 

Formal limit, if it exists, is unique. Formal convergence is stronger than hybrid: 

c = flim„ c(n) => [c] = hlim„ c(n) 

and we presented sequence with hybrid limit but without formal limit. Further, 
flim„ c\n = c and hlim„ c\n = [c] for every c G D. 

Proposition 2.8 If c(n) and d{n) are sequences of terminating decimals that 
have hybrid limits and c(n) < d(n) for every large n, then 

hlim„ c(n) < hlim„ d(n). 

Proof. We denote [c] = hlim n c(n), [d] — hlim„ d{n) and assume for contrary 
that [c] > [d]. Thus c > d and, since c ^ d, there exist k,N G N such that 
c\n — d\n > 10~ fc for n > N. As c(n) — c\n — > and d{n) — d\n — > 0, we see that 
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there is an AT' such that c(n) — d(n) > 10 k 1 for n > N'. This contradicts the 
assumption on c(n) and d(n). □ 



Recall that a sequence d(n) of terminating decimals is Cauchy if for every 
k G N there is an N G N such that |d(m) — < 10 _fc whenever m,n > N. 

Proposition 2.9 Let d(n) be terminating decimals. 

1. If the sequence d(n) is monotone and bounded, then it has a formal limit. 

2. If the sequence d(n) is Cauchy, then it has a hybrid limit. 

Proof. 1. We suppose that d(n) is nondecreasing, the other case is similar. 
Using Proposition ^. 4[ we set d — sup ra d(n). We prove that fhm„ d(n) = d. If d 
is the larger element of a jump (i.e., d > 0, d G To or d < 0, d G Tg), then in fact 
d = max,, d(n) and d(n) from an index on equals d. Else for any given I G Z 
there is a decimal e such that e < d and e; = <i, for z > I. Then e < d(N) < d for 
some N and it follows that d(N)i — e, = rf; for i > L Hence, by monotonicity, 
the same is true for every d(n) with n > N. Thus d(n) formally converge to d. 

2. The proof is the same as in the classic case for metric limit and real 
numbers. Sequence d(n) is Cauchy and therefore bounded. It contains a mono- 
tone subsequence d(rik) (every infinite sequence in a linear order has an infinite 
monotone subsequence). By part 1, ffinifc d(rik) = c for some c e D. Thus 
hlinifc d(nk) = [c]. By the Cauchy property of d(n) and hlim^ d{nk) — [c], for 
given I 6 N there is an N 6 N such that m > n > N implies \d(m) — d(n)\ < I0~ l 
and k > N implies \d(rik) — c\k\ < 10 - '. Thus for every k > N we have that 
\d{k) - c\k\ < \d{k) - d(n k )\ + \d{n k ) - c\k\ < 2 • 10~ fc . Thus d(k) - c\k -> and 
hlim fe d(k) =[c]. □ 



Proposition 2.10 The sum a + (3 and product a(3 is defined for every pair of 
real numbers a, (3 € R. 

Proof. If c(n) and d(n) are two Cauchy sequences of terminating decimals, then 
c(n) + d(n) and c(n) • d(n) are Cauchy as well (to see the latter use the above 
rearrangement). For any decimals c 6 a and d £ j3 the sequences c|n and d|n 
are Cauchy, in fact 

| c|m — c|n | < 10 _n for m > n 

and similarly for d\n. Thus c|n + d|n and c\n ■ d\n are Cauchy, have hybrid 
limits by 2 of Proposition 12 . 91 and define the sum a + j3 = [c] + [d] and product 
a-/3=[c]- [d]. ' ' □ 

Note that for terminating decimals c, d one has [c] + [d] = [c+d] and [c][d] = [cd]. 
For every fc G Z and d G D one has [I0 k ][d] = [10 k d]. 

We mentioned in the Introduction that it is necessary to justify exchanges 
of limits with sums and products; for hybrid limits it is easy. 
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Proposition 2.11 Let c(n) and d(n) be sequences of terminating decimals that 
have hybrid limits. Then so have the sequences c(n) + d(n) and c(n) ■ d{n) and 



hlim„ (c(n) + d(n)) 
hlim„ (c(n) • d(n)) 



hlim„ c(n) + hlim„ d{n) 
hlim„ c(n) • hlim„ d(n). 



Proof. Let hlim„ c(n) = [c], hlim n d(n) = [d] and [c] + [d] = [e]. Then 
c(n) — c|n — > and d(n) — d\n — > 0, which sums to c(n) + (i(n) — (c|n + d|n) — > 0. 
By the definition of addition, c\n + d\n — e\n — > 0. Summing again we get 
c{n) + d(ri) — e\n — > 0. Thus hlim„ (c(n) + d{n)) — [e]. 

As for the product, let [c] • [d] = [e]. We noted above that c(n) — c\n — > 
and d(n) — d\n — > implies c(n) • d(n) — c|n • d\n — > 0. By the definition of 
multiplication, c\n ■ d\n — e\n — > 0. Summing we get c(n) • d(n) — e\n — ► 0. Thus 
hlim„ (c(n) • d(n)) = [e]. □ 

Conclusion of the proof of Theorem 12.21 We need to prove that addi- 
tion and multiplication on M. are commutative and associative operations with 
distinct neutral elements and 1, have inverse elements (except for multipli- 
cation), satisfy the distributive law, and that (M, <) is a complete linear order 
to which addition of an element and multiplication by a positive element are 
increasing functions. 

We first prove that (R, +, •) is an integral domain. It is clear that [0] = {0} 
and [1] = {0.999 . . . , 1} are neutral elements to -I- and • and that the commu- 
tativity of + and • is transferred by hybrid limits from To to R. Also, it is 
clear that every [c] G R has [— c] as its additive inverse. However, associativity 
and distributivity need Proposition 12.111 For example, to prove that addition 
is associative, we take [c], [d] and [e] in R and, using Propositions 12.51 and 12.111 
get that indeed 

([c] + [d]) + [e] = hlim n (c\n + d\n) + hlim n e\n (definition of addition) 



In the same way it follows by Propositions 12.51 and 12.111 that multiplication is 
associative and that the distributive law holds. If [c][d] = [0], then c\n-d\n — > 0, 
which implies (as both c\n and d\n is monotone) that c\n — > or d\n — > 0, that 
is, [c] = [0] or [d] = [0]. Thus (R, + , •) is an integral domain. 

We consider the properties of <. In Proposition 12.41 we have already proven 
that (R, <) is a complete linear order. We need to show that 



We have c < d and therefore c\n < d\n for every n and thus (by Proposition ^. 5| 
c|ri + e|n < d\n + e|n for every n and any real number [e]. Taking hybrid 



= h\im n ((c\n + d\n) + e\n) ( Proposition 12 . 1 lj) 
= hlimn (c|n + (d\n + e\n)) (Proposition 12. 5p 



= hlim„ c\n + hlim n (d\n + e\n) (Proposition 12. llj) 
= [c] + ([d] + [e]) (definition of addition). 



[c] < [d] => [c] + [e] < [d] + [e] and, for [e] > 0, [c][e] < [d][e]. 
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limits, by Proposition 12.81 we get that [c] + [e] < [d] + [e]. But the equality 
[c] + [e] = [d] + [e] is impossible as cancelling [e] (i.e., adding [— e] to both sides 
and applying associativity of +) would give [c] = [d]. Hence [c] + [e] < [d] + [e]. 
For multiplication we argue similarly; the cancellation of [e] in [c] [e] = [d] [e] 
follows from the fact that (R, +, •) is an integral domain. Thus (R, + ,- ,<) is an 
ordered ring. 

It remains to show that multiplication has inverses. 
Proposition 2.12 For every nonzero decimal c there is a decimal d such that 

[c] ■[d] = [l} = {0.999... ,1}. 

Proof. It suffices to construct a Cauchy sequence of terminating decimals d(n) 
with the property that c\n ■ d(n) — 1 — > 0. Setting [d] = hlim„ d(n) and taking 
hybrid limits, we get by Proposition 1 2 . 1 1 1 that 

[0] — hlim„ (c|n • d{n) — 1) = hlim„ c\n ■ hlim„ d(n) — hlim„ 1 = [c] • [d] — [1] 

and (by the already established properties of +) [c] • [d] = [1]. 

We may assume that c is positive. For given n £ N, we consider decimals 
c' = 10™ • c\n and 10^ in Z , where N e N is selected so that lO^ - " > d . 
Dividing in Z , we have 10^ = c'd + e with d, e £ Z , d > and < e < c'. 
Multiplying by lO - ^, we get in T the equality 

1 = {10- n c')(lQ- N+n d) + l0~ N e 
= c\n ■ {10- N+n d) + 1Q- N e 
—: c\n ■ d(n) + e{n). 

As 

< e[n) = 1Q- N e < lO^c' < lO - ", 

we have 1 — c\n ■ d(n) — e(n) — > 0. We show that d(n) is Cauchy. Because c > 0, 
there exist fceZ and N € N such that 

c|n > I0 k for n> N. 

From < c\n ■ d(n) = 1 — e(n) < 1 we get 

< d(n) < 10- k for n > N. 

We take m > n > N and write c\m — c\n + S and d(m) = d(n) + A. Recall that 

< 8 < 10" n . 

Subtracting 1 = c|n ■ d(n) + e(n) from 1 = (c\n + S)(d(n) + A) + e(m) and 
rearranging, we get 

10 fe |A| < \c\m- A | = \e(n) - e{m) - 8 ■ d{n)\ 

< \e(n)\+\e(m)\+\S-d(n)\ 

< 2 • 10"" + 10"™" fc . 
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Thus |A| = \d(m) - d(n)\ < 2 • W- n - k + 1CT™- 2A; for m > n > N and we see 
that d(n) is Cauchy. □ 

This concludes the proof of Theorem 12.21 

The proof of Theorem 12.31 We start with an arithmetic characterization 
of ultimately periodic real numbers. For a, 6 G N, we let 9^0^ denote the 
decimal 99 ... 900 ... G Z with a 9's and b 0's and similarly for 9< a ). 

Proposition 2.13 A real number a 6i is in [P] if and only if [9^0^} a is in 
[Zo] for some a G N. 

Proof. Suppose that a — [c] for c G P with period p. Then there is an m G Z 
such that for every k G N we have [9 (fep) ][c] = [10 fcp ][c] - [c] = [10 fcp c] - [c] = [d] 
where di = for i < m. Taking A: large enough so that a = kp > |m|, we get 
[10 a ][9( a )][c] = [9( Q )0( a )][c] in [Zo]. 

Now suppose that for a G N and c G D we have [9( a W a >][c] G [Z ]; we may 
assume that c is positive. It follows that there is a k G N and a terminating 
decimal e such that [10 a c] — [c] = [e] and a — for i < —k. For I > k we write 
(10 a c)|Z = (10 a c)|fc + 5(0 and c|Z = c|fc + A(l). It is clear that 6(1), A(l) G T 
and are nonnegative, smaller than 10 _fc , and, since 10 a c ^ To or c G" Tg, there 
is an m > k such that for every I > m is 6(1) > 10~ m or for every I > m is 
A(l) < 10- k - 10~ m . Thus for every / > m, 

| (I0 a c)\l - c\l - ((10 a c)|fc - c\k) | = \6(l) - A(/)| < 10- fc - 10- ro . 

If (10 Q c)|fc - c\k ^ e, then | (10 a c)|fc — cjfc — e | > 10" fc and, using that e - 
((10 a c)|Z — c\l) — > 0, by summing we get for large Z a contradiction with the 
displayed inequality. Hence (10 a c)\k — c\k = e. Since we can increase fc, we get 
(10 a c)|fc' — c\k' = e for any k' > k. So (10°c)j = Cj for every i < —k, which 
means that c G P with period a. □ 



Lemma 2.14 For every nonzero number r in Zo there is an a G N such that 
g( a )()( a ) is divisible by r. 

Proof. We split r — st where s has only prime factors 2 and 5 and t is coprime 
with 10 = 2 • 5. By Euler's generalization of the little theorem of Fermat, for 
p = ip(\t\) is 10 fep — 1 = 9^ fcp ' divisible by t for any k G N. Taking fc large enough 
so that I0 kp is divisible by s, we get that lO^lO^ - 1) = 9( fep W fcp ) is divisible 
by st = r. □ 

We show that [P] is the prime field of R. The prime field is formed by the 
elements a/b, where 6^0 and a and b are of the form ±([1] + [1] + . . . + [1]). 
It is clear that these are exactly the elements in [Zo]. If a G [P], then by 
Proposition [233] [9( Q W a ']a = for some a G N and j3 G [Z ]. Thus 

a = /V[9 (Q) (a) ] 
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is in the prime field. So [P] is contained in the prime field. 

Now it suffices to show that [P] is a subfield of R. Let a, (3 £ [P] be any 
elements. Using Proposition 12.131 and Lemma 12.141 we take numbers a, b, c in 
N so that [9^0^}a e [Z ], [9< 6 W 6 >]/3 e [Z ] and [9< c )o( c )] is divisible by the 
product [9 (a )0(°)][9W0 (h) ]: [9^0^] = [9^0^][9^0^]r for some r £ [Z ]. 
Then 

[9 (c) (c) ](a ±P) = (r[9 (b) (b) ])([9 (a) (a) ]a) ± (r[9 (a) (a) ])([9 (b) (b) ]/3) 
is in [Zq] and so is 

[9 (c) (c) ](a/?) = r([9 (a) (a) ]a)([9 (b) (b) ]/3). 

Thus, by Proposition 12. 131 a ± (3 and a/3 are in [P]. If a ^ 0, we take a d £ N 
so that [9( d W d >] is in [Z ] divisible by [9( a W a >]a and get that 

[9^0 (d) ](l/a) = [9 (d) (d) ][9 (a) (a) ]/[9 (a) (Q) ]a 

is in [Zq], which means that 1/a £ [P]. Thus [P] is a subfield of R and must be 
equal to the prime field. This concludes the proof of Theorem 12.31 

Addition and multiplication of decimals. We showed how to add and 

multiply real numbers [c] and [d] without worrying what is c + d and cd. For 
completeness, we prove that formal limits define addition and multiplication 
on D correctly. We already established this with the exception of subtraction. 
Indeed, for any c, d £ D the sequence c|n • d\n is monotone and bounded and 
so flim n (c|n • d\n) exists by part 1 of Proposition 12.91 and defines correctly the 
product cd. This works for c + d defined by flim„(c|n + d\n) as well unless c 
and d have opposite signs. If this happens, c|n + d\n may not be eventually 
monotone and we need a finer criterion of formal convergence. 

Proposition 2.15 Let d(n) be a sequence of terminating decimals. Suppose 
that for every k £ N there exist an N £ N and a pair c(k) < c(k)' of terminating 
decimals such that c{k)i = c(k)' i = for i < —k, c(k)' — c(k) = 10~ fc and c(k) < 
d(n) < c{k)' for every n > N. The sequence d(n) then formally converges. 

Proof. Observe that 

c(l) < c(2) < ...<c<d < ...< c(2)' < c(l)' 

where c = sup fc c(fc) and c' = inf/j. c(k)' . By part 1 in Proposition ^. 91 

flinifc c(fc) = c and flim/j c{k)' = c . 

Suppose first that c = c! . For given I £ N there is a A; such that c{k)i = c{k)\ = 
c, = c[ for every i > —I. For this k we have c(k) < d(n) < c(k)' for every 
n > N. It follows that d(n)i = Ci = c\ for every i > —I and n > N. Thus 
flim„ d(n) = c = d ' . 



14 



The remaining case is that c < d but c ~ d ' ■ We assume that c and d are 
negative and c E To, d <E Tg, the other case is similar. It follows that c(k) is 
eventually constant, equal to c. For given I S N we fix a k such that c(k) = c 
and c(k)' i = for every i > — Z. For this /c we have c < d(n) < c(k)', hence 
d < d(n) < c(k)' , for every n > N. It follows that d(n)i = d i for every i > —I 
and n> N. Thus nim„ d(n) = c'. □ 

Note that the criterion becomes invalid when the condition c(k) < d(n) < c(k)' 
is weakened to c(k) < d(n) < c(k)' . 

Proposition 2.16 Let c, d £ D oe decimals. The formal limits 

flim n (c|n + and fhm„ (c|n • d|n) 
exist and define the sum and product of decimals c + d and c ■ d. 

Proof. As we mentioned, it only remains to deal with subtraction c — d of 
positive decimals c and d. We may assume that q ^ dj for infinitely many i 
because else c|n — d|n is eventually constant and ffim„ (c|n — d|n) exists trivially. 
Let k G N be given. For I > k we write c\l = c|fe + (5(0 and d|Z = d\k + A(l). 
The terminating decimals 5(1) and A(Z) are nonnegative and smaller than 10^. 
By the assumption on c and d, there is an I' > k such that 5(1') — A(l') ^ 0. 
We assume that 5(1') — A(Z') < 0, the other case is similar. It follows that 
5(1') - A(Z') < -10"'' and 5(1) - A(Z) < for every / > I'. Setting c(k) = 
c|fc - d|fc - lCT fe , c(fc)' = c|fc - d\k and AT = I', for n > TV we have 

c(fc) < c\n - d\n = c\k - d\k + 5(n) - A(n) < c(k)' . 

Thus ffim„ (c\n — d\n) exists by Proposition ^. 151 □ 

By properties of formal and hybrid limits we have equalities 

[c] + [d] = [c + d] and [c][d] = [cd] 

for every pair of decimals c, d G D, not just for terminating ones. Addition 
and multiplication of decimals is trivially commutative and by Theorem 12.21 
also associative and distributive up to equivalence, that is to say, (c + d) + e ~ 
c + (d + e), (cd)e ~ c(de) and c(d + e) ~ cd + ce for every c, d, e S B. In fact, 
for multiplication it holds more strongly that 

(cd)e = c(de) for every c, d, e € D. 

This can be established by a version of Proposition 12.111 for formal limits — 
details are left for the interested reader. 

References and remarks. The oldest model of real numbers uses cuts on the 
set of rational numbers and is due to Dedekind who conceived it, by his own 
words, on November 24, 1858 (Dedekind [12] where he after 14 years outlines his 
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approach). Modern version is, e.g., in Pugh [331 Chapter 1.2] or Rudin [351 P- 
17]. The first rigorous theory of irrational numbers to appear in print was the 
construction based on (equivalence classes of) Cauchy sequences of rationals, 
due to Meray [3D] in 1869 and later Cantor [pj and Heine [25]. For modern 
treatment see, e.g., Tao [33J Chapter 5]. 

Another model of K developed Weierstrass in his lectures (starting 1865/66). 
He published nothing on his theory but it was disseminated through his students, 
some of which (Biermann, Hettner, Hurwitz, Killing, Kossak, Pasch, Pincherle — 
see [35J and [3PJ P- 46]) produced written accounts on his approach; thorough 
treatment (based on Weierstrass' courses in 1872 and 1884) is due to Dantscher 
[TT] . Similar theory of irrational numbers was developed by Stolz [35], [4"0] , 
This approach is often dubbed as Weierstrass or Weierstrass-Stolz model of R 
based on decimal expansions (Gamelin [18) . Kudryavtsev [28] and elsewhere). 
We looked in a few primary sources (Stolz [39] and Stolz and Gmeinder [40] . 
Dantscher [TT] was available to us only in the detailed contemporary review of 
Miller [3T], see also Snow [35]) and, unsurprisingly, did not find there the modern 
decimal model of R; decimal expansion is mostly treated not as one object, 
actually infinite string of symbols or numbers, but via partial sums as a sequence 
of rational approximations, which is not very different from the Meray-Cantor- 
Heine model. But in at least one case some features of decimal arithmetics 
appear — Stolz and Gmeinder [J0J pp. 48-50] give, in effect, an algorithm for 
calculating the sum c + d of two infinite decimals, justify subtraction and prove 
the associativity (up to equivalence) of addition of infinite decimals (in the form 
that -3/10™ < {{a + b) + c}-{a+{b + c)} < 3/10" foreveryn). However, after 
that they write: "It would come out even more tedious to explain multiplication 
of two real numbers. For building the four arithmetic operations one prefers now 
to switch to the Cantor's or Weierstrass' theory of irrational numbers." 

We found several outlines or sketches of the decimal model of R, namely 
Courant and Robbins [10], Gamelin [18], Gowers [2pJ Chapter 4], [21], and [22], 
Kudryavtsev [25J , and Richman [35J but no really detailed and rigorous account. 
This was a motivation to present our detailed decimal model of R in Section 2. 

There are other constructions of real numbers in the literature: Eudoxus 
real numbers (Arthan [3J, A'Campo [5J, Douglas et al. [T3J and Street [UJ), 
the construction of de Bruijn ([6]), factorization of a set of power series (Knuth 
and Pratt [27]), an algebraic approach of Faltin ct al. (15 ). For more on real 
numbers see the books of Ebbinghaus et al. [14] and especially (unfortunately 
only for those reading Slovak) Bukovsky [7]. 

3 Irrationality of \/2 in the decimal model of R 

We use the decimal models of R and Q as developed in the previous section 
in Theorems 12.21 and 12.31 and their proofs. In particular, recall that for every 
c, d G To, [c] = [d] implies c = d and that 

[c] + [d] = [c + d] and [c][d] = [cd]. 
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The displayed equalities hold in fact even for every c, d £ D but this generality 
will be needed only in the proof of Proposition 13.31 and not in Proposition 13.11 
Proposition 13.11 is on the arithmetics in R = [D] based on hybrid limits and 
Proposition 13.31 is on the arithmetics in D based on formal limits. 

Proposition 3.1 In the decimal model of real numbers, the equation x 2 = [2] 
has no solution in the set of ultimately periodic real numbers [P]. 

Proof. We assume for contrary that some a £ [P] satisfies a 2 — [2] . By Propo- 
sition we have fi^O^a = [e] for some a £ N and e £ Z . Thus 

[e] 2 = ([9( a ¥ a) ]a) 2 = [9 (a) (a) ] 2 [2] and e 2 = (9 (a) (a) ) 2 2. 

As we know from Lindstrom's proof, the last equality is impossible because the 
last nonzero digit on the left is 1,4,5,6 or 9 but the last nonzero digit on the 
right is 2. □ 

There is nothing special about 2, the same argument works if the last nonzero 
digit is a quadratic non-residue modulo 10. We obtain the following result. 

Proposition 3.2 If c is a positive terminating decimal whose last nonzero digit 
is 2, 3, 7 or 8 then the equation x 2 = [c] has no solution in the set of ultimately 
periodic real numbers [P]. 

But we want to get the stronger result that x 2 = 1.999 . . . has no solution 
in P. The author once devised the following "proof". If d 2 = 1.999 . . . for some 
d £ P with period p, then 10 p d — d = e £ To and 

e 2 = (I0 p d - d) 2 = d 2 (10 p - l) 2 = 1.999 . . . • (9 (p) ) 2 . 

It is easy to see that any product ab with a £ T 9 and b £ T \{0} lies in Tg 
and therefore the right side shows that e 2 lies in Tg. This is clearly impossible 
because e £ T , e 2 £ T and T and Tg are disjoint. Hence x 2 = 1.999 . . . has 
in D no ultimately periodic solution. 

Unfortunately, this quick argument is fallacious because of the failure of the 
distributive law in D. The second equality sign in the calculation is correctly 
equivalence, (10 p d — d) 2 ~ <i 2 (10 p — l) 2 , and the calculation shows only that 
e 2 ~ 1.999... -(90)) 2 , which is no (immediate) contradiction. Indeed, the same 
argument would prove that x 2 = 3.999 . . . has no solution in P either, which is 
in conflict with the existence of the ultimately periodic solution x = 1.999 .... 

Proposition 3.3 The equation x 2 = 2 has no solution in the set of decimals 
D and the equation x 2 = 1.999 . . . has no solution in the subset of ultimately 
periodic decimals P. 

Proof. We assume the contrary that some decimal d satisfies d 2 = 2. Since 
flim„ (d\n) 2 = 2, there is an N £ N such that for n > N we have ((d\n) 2 )i = 
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for i > and ((d\n) 2 ) = 2. We fix an m larger than N. Then ((d\m) 2 )j > for 
some j < as 2 is not a square of any terminating decimal (because of the last 
nonzero digits). Since the sequence (d\n) 2 is nondecreasing, for every n > m 
we have (<i|n) 2 > (d\m) 2 and therefore {{d\n) 2 )i > for some i, j < i < 0. But 
flim n (d\n) 2 — 2 implies that {{d\n) 2 )i = for every i, j < i < 0, and large n, 
which is a contradiction. 

The unsolvability of the equation x 2 = 1.999 ... in P follows immediately 
from Proposition 13. II if some d € P satisfies d 2 — 1.999 . . ., then [d] 2 = [d 2 ] = 
[1.999 ■ ■ •] = [2] in R, which was shown to be impossible. □ 

The completeness of (D, <) provides the aperiodic solution(s) of x 2 = 1.999 . . . , 

x = ±1.414213562373095048801688724209698 . . . 
For algorithms calculating this decimal see Sebah and Gourdon [37] , 
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